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Abstract 

We scan the parameter space of the NMSSM for the observability of a Higgs boson 
at the LHC with 300 fb _1 integrated luminosity per detector, taking the present LEP 
constraints into account. We focus on the regions of parameter space for which none 
of the usually considered LHC detection modes are viable due to the fact that the only 
light non-singlet (and, therefore, potentially visible) Higgs boson, h, decays mainly to 
two CP-odd light Higgs bosons, h — > aa. We simulate the WW — > h — > aa detection 
mode. We find that this signal may be detectable at the LHC as a signal/background ~ 
(500—1000) / 300 bump in the tail of a rapidly falling mass distribution. If further study 
gives us confidence that the shape of the background tail is predictable, then we can 
conclude that NMSSM Higgs detection at the LHC will be possible throughout all of 
parameter space by combining this signal with the usual detection modes previously 
simulated by ATLAS and CMS. We also show that this WW — > h — » aa signal will 
be highly visible at the LC due to its cleaner environment and high luminosity. We 
present a study of the production modes and decay channels of interest at the LC. 
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1 Introduction 



One of the most attractive supersymmetric models is the Next to Minimal Supersymmet- 
ric Standard Model (NMSSM) dEHHm, which extends the MSSM by the introduction 
of just one single superfield, S. When the scalar component of S acquires a TeV scale 
vacuum expectation value (a very natural result in the context of the model), the super- 
potential term SH u H d generates an effective fiH u H d interaction for the Higgs doublet su- 
perfields. Such a term is essential for acceptable phenomenology No other SUSY model 
generates this crucial component of the superpotential in as natural a fashion. Thus, the 
phenomenological implications of the NMSSM at future accelerators should be consid- 
ered very seriously One aspect of this is the fact that the h, if, A, H ± Higgs sector of the 
MSSM is extended so that there are three CP-even Higgs bosons (/ii,2,3/ mhi < mh 2 < mh 3 ), 
two CP-odd Higgs bosons (ai,2, rn ai < m a2 ) (we assume that CP is not violated in the 
Higgs sector) and a charged Higgs pair (h^). An important question is then the extent 
to which the no-lose theorem for MSSM Higgs boson discovery at the LHC (after LEP 
constraints) is retained when going to the NMSSM; i.e. is the LHC guaranteed to find at 
least one of the ^1,2,3/ «i,2/ h~^l 

One of the key ingredients in the no-lose theorem for MSSM Higgs boson discovery 
is the fact that decays of the SM-like Higgs boson to AA are only possible if m A is quite 
small. This is due to the fact that the masses of the h, H, A, H ± of the MSSM are closely 
tied to one another. In particular, can only be small if m h is also relatively small 
and m H ± ~ m w . In this region, it is the H that is SM-like. The small-m^ region is either 
excluded by LEP because of non-observation of hA production, or will be detectable at the 
LHC via observation of t — ► bH ± decays in tt production. Most of the still unconstrained 
portion of MSSM parameter space corresponds to the decoupling region of ~ m H ~ 
m H ± > 120 - 130 GeV in which the h is SM-like and m h <120 - 130 GeV so that h —> AA 
decays are forbidden. In the NMSSM, these strong mass relations are weakened. The 
hi or hi can be SM-like without the a\ necessarily being heavy. As a result, h\ — > a±ai 
or hi — > a\d\ decays can be prominent and fall outside the scope of the usual detection 
modes for the SM-like MSSM h on which the MSSM no-lose LHC theorem largely relies. 

In Ref . 0, a partial no-lose theorem for NMSSM Higgs boson discovery at the LHC 
was established. In particular, it was shown that the LHC would be able to detect at 
least one of the NMSSM Higgs bosons (typically, one of the lighter CP-even Higgs states) 
throughout the full parameter space of the model, excluding only those parameter choices 
for which there is sensitivity to the model-dependent decays of Higgs bosons to other 
Higgs bosons and/or superparticles. Here, we will retain the assumption of a heavy 
superparticle spectrum and address the question of whether or not this no-lose theorem 
can be extended to those regions of NMSSM parameter space for which Higgs bosons can 
decay to other Higgs bosons. We find that the parameter choices such that the "standard" 
discovery modes fail would allow Higgs boson discovery if detection of h — > aa decays is 
possible. (When used generically, the symbol h will now refer to h = hi, h 2 or h 3 and the 
symbol a will refer to a = ai or a 2 ). Detection of h — > aa will be difficult since each a will 
decay primarily to either bb (or 2 jets if m a < 2mb) and r + r~, yielding final states that will 
typically have large backgrounds at the LHC. Further, a light a can only be detected on 
its own at the LHC if it has a highly enhanced Yukawa coupling to b (or t) quarks. 

Let us begin by reviewing our earlier procedure. For each point in parameter space 
not ruled out by constraints from LEP on Higgs boson production, e + e — > Zh Il6ll or 
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associated production e + e~ — > /ia Q, we calculate the Higgs boson masses and decay 
branching ratios including all the relevant radiative corrections. We eliminate parameter 
choices for which one Higgs boson can decay to two other Higgs bosons or a vector boson 
plus a Higgs boson. For the remaining regions of parameter space, we then estimate the 
statistical significances (computed as N SD = S/ \fB for a given mode) for all Higgs boson 
detection modes so far studied at the LHC HSll^ HTOllTTI . These are (with I = e, /i) 
t)gg -> h/a -> 77; 

2) associated Wh/a or tth/a production with r yy£ ± in the final state; 

3) associated tth/a production with h/a — > bb; 

4) associated bb h/a production with h/a — > t + t~; 

5) gg^h^ ZZ& -> 4 leptons; 

6) gg^h^ WWW -> ^"i/P; 

7) WW -> /i -> r+r-; 

8) WW -> /i -> WWW. 

The outcome is that, for an integrated luminosity of 300 fb 1 at the LHC, there are still 
regions in the parameter space with < 5a expected statistical significance for all Higgs 
boson detection modes so far studied in detail by ATLAS and CMS, i.e. modes 1) - 6) but 
not the W-fusion modes 7), 8). On the other hand, the expected statistical significance for 
at least one of the non-W-fusion modes is always above 3.6a at 300 fb~ , and the statistical 
significance obtained by combining (using the naive Gaussian procedure) all the non-W- 
fusion modes is at least 4.8a. However, we found that all such cases are quite observable 
(at > 10.1a) in one of the W-fusion modes (using theoretically estimated statistical sig- 
nificances for these modes). For all points in the scanned parameter space, statistical 
significances obtained by combining all modes, including W-fusion modes, are always 
>10.7a. Thus, NMSSM Higgs boson discovery by just one detector with L = 300 fbT 1 
is essentially guaranteed for those portions of parameter space for which Higgs boson 
decays to other Higgs bosons or supersymmetric particles are kinematically forbidden. 

In this work, we investigate the complementary part of the parameter space, where 
at least one Higgs boson decays to other Higgs bosons. To be more precise, we require at 
least one of the following decay modes to be kinematically allowed: 

i) h — > h'h! , ii) h — > aa , Hi) h — *> h , iv) h — > aZ , 
v ) h — > h^W T , vi) a' — > ha , vii) a — » hZ , viii) a — > h ± W T . (1) 

The observability at the LHC of Higgs bosons that decay in these ways has not been stud- 
ied, except for some particular MSSM cases that are not applicable for the masses and 
Yukawa coupling strengths that are most relevant in the NMSSM context. After search- 
ing those regions of parameter space for which one or more of the decays i) — viii) is 
allowed, we found that the only subregions that would not be excluded by modes 1) - 8) 
correspond to NMSSM parameter choices for which (a) there is a light CP-even Higgs bo- 
son with substantial doublet content that decays mainly to two still lighter CP-odd Higgs 
states, h — > aa, and (b) all the other Higgs states are either dominantly singlet-like, im- 
plying highly suppressed production rates, or relatively heavy, decaying to tt or to one of 
the "difficult" modes i) — viii). In such cases, it seems evident that the best opportunity 
for detecting at least one of the NMSSM Higgs bosons is to employ WW — > h production 
and develop techniques for extracting a visible signal for the h — > aa final state. We have 
performed a detailed simulation of this signal and find that its detection may be possible 
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after accumulating 300 fb _ in both the ATLAS and CMS detectors. The nominal statis- 
tical significance achieved is > 30cr with S/B > 1. However, the signal only emerges on 
the tail of a rapidly falling background. A thorough understanding of and confidence in 
the shape of this background tail would be needed to claim detection of the h — > aa sig- 
nal. If further study and, ultimately, actual measurements of the background outside the 
signal region result in such confidence then the combined search for this signal and the 
signals for Higgs bosons in the previously studied modes 1) - 8) would result in detec- 
tion of at least one of the NMSSM Higgs bosons throughout the entire NMSSM parameter 
space (assuming no large branching ratios for Higgs boson decays to sparticles). Having 
reached this conclusion, we turn to a study of the role that a future LC will play in explor- 
ing the NMSSM Higgs sector for those parameter choices for which detecting the h — > aa 
signal is so critical. We find, regardless of what happens at the LHC, that the LC is guar- 
anteed to find one or more Higgs bosons, including detection of the h — > aa signal in 
several different production channels, leading to the possibility of measuring the h — ► aa 
branching ratio. 

In section |2j we review the basic set up of the NMSSM as well as our scanning pro- 
cedure for the NMSSM parameter space. In section |3J we give a selection of benchmark 
points for which all statistical significances at the LHC are below 5a after employing the 
standard discovery modes 1) - 8). We then discuss possible new discovery channels in 
section |U together with our simulation results both at the LHC and the LC. Conclusions 
are given in section |5j 



2 NMSSM Parameters 

We consider the simplest version of the NMSSM (1113111/ where the term fiHiH 2 in the 
superpotential of the MSSM is replaced by (we use the notation A for the superfield and 
A for its scalar component field) 

\H X H 2 S + ^S 3 , (2) 

so that the superpotential is scale invariant. We make no assumption on "universal" soft 
terms. Hence, the five soft supersymmetry breaking terms 

m\H\ + m\ 2 H\ + m 2 s S 2 + \A X H X H 2 S + ^A K S 3 (3) 

are considered as independent. The masses and/ or couplings of sparticles are assumed 
to be such that their contributions to the loop diagrams inducing Higgs boson production 
by gluon fusion and Higgs boson decay into 77 are negligible. In the stop sector, which 
appears in the radiative corrections to the Higgs potential, we chose the soft masses rug = 
rriT = M susy = 1 TeV, and vary the stop mixing parameter 

As in the MSSM, the value X t = - so called maximal mixing - maximizes the radiative 
corrections to the Higgs boson masses, and we found that it leads to the most challenging 
points in the parameter space of the NMSSM. 
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Assuming that the Higgs sector is CP conserving, the independent parameters of the 
model are thus: A, k, vr? H , , m 2 Sl A x and A K . For purposes of scanning and analysis, it 
is more convenient to eliminate m 2 H , rn 2 H2 and m| in favor of M z , tan/? and fi c g = X(S) 
through the three minimization equations of the Higgs potential (including the dominant 
1- and 2-loop corrections 1121 ) and to scan over the six independent parameters 

A, k, tan P, /x e ff, A\, A K . (5) 

We adopt the convention A, k > 0, in which tan (3 can have either sign. The absence 
of Landau singularities for A and n below the GUT scale (~2x 10 16 GeV) imposes upper 
bounds on these couplings at the weak scale, which depend on the value of h t and hence 
of tan j3 |HJ|3|. Using mf ole = 175 GeV, one finds A max ~ 0.69 and K max ~ 0.62 for interme- 
diate values of tan (3. In addition, we require |/i e ff| > 100 GeV; otherwise a light chargino 
would have been detected at LEP. (The precise lower bound on \fM e s\ depends somewhat 
on tan (3 and the exact experimental lower bound on the chargino mass; however, our 
subsequent results do not depend on the precise choice of the lower bound on \/i e s\-) 

We have performed a numerical scanning over the free parameters, which were ran- 
domly chosen in the following intervals: 

0.1 < A < A max , 0.001 < k < /« max , 2.5 < | tan/?| < 10 , 

100 GeV < |/i eff | < 500 GeV , < \A X \ < 500 GeV , < \A K \ < 500 GeV . (6) 

For each point, we computed the CP-even and CP-odd Higgs boson masses and mixings, 
taking into account radiative corrections up to the dominant two loop terms, as described 
in 1121 . The five mass eigenstates are denoted by h{ (i = 1, 2, 3) and aj (j = 1, 2) for CP- 
even and CP-odd eigenstates respectively, with masses m h . and m a . in increasing order. 
We have eliminated parameter choices that lead to violation of the LEP constraints on 
Higgs boson production e + e — > Zhi (i = 1, 2 or 3) [6]. These constraints give an upper 
bound on the ZZhi reduced coupling, R ir as a function of . (The reduced coupling Ri is 
defined as the coupling ZZhi divided by the corresponding Standard Model coupling and 
is the equivalent of sin(/3 — a) in the MSSM.) We have also kept only parameters leading 
to consistency with the exclusion limit from LEP on Higgs bosons associated production 
e + e~ — >• hi<ij [7J. This provides an upper bound on the ZhiCLj reduced coupling, R' ijr as a 
function of + m aj for — m aj Q- (The Zhidj reduced coupling is the equivalent of 
cos(P — a) in the MSSM.) Finally, we calculated the charged Higgs boson mass m h ± and 
required m h ± > 155 GeV, so that t — > h^b would not be seen. 

In order to probe the complementary part of the parameter space as compared to the 
scanning of Ref. @, we required that at least one of the "difficult" decay modes i) — viii) 
displayed in the previous section is allowed. For each Higgs state, we have calculated the 
branching ratios in the usual decay modes as well as in the difficult modes i) — viii) when 
kinematically allowed. In doing so, we have included all relevant higher-order QCD 
corrections |13| using an adapted version of the FORTRAN code HDECAY 1141 . We then 
estimated the expected statistical significances at the LHC in all Higgs boson detection 
modes 1) - 8) as described in the next section. 
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3 Invisible Points at the LHC 



From the known couplings of the NMSSM Higgs scalars to gauge bosons and fermions it 
is straightforward to compute their production rates in gluon-gluon fusion and various 
associated production processes, as well as their partial widths into 77, gauge bosons 
and fermions, either relative to a SM Higgs scalar or relative to the MSSM h, H and/ or A. 
This allows us to apply "NMSSM corrections" to the processes 1) - 8) of section[TJ These 
NMSSM corrections are computed in terms of the following ratios: For the scalar Higgs 
bosons, Ri is the ratio of the coupling of the hi to vector bosons as compared to that of a 
SM Higgs boson, and ti, h are the corresponding ratios of the couplings to top and bottom 
quarks. (The coupling ratio for r leptons is always the same as for b quarks, i.e. Ti = bi). 
Note that we always have \Ri\ < 1, but ti and bi can be larger, smaller or even differ in sign 
with respect to the SM. For the CP-odd Higgs bosons, there is no tree-level coupling to the 
VV states; t'j and bj are defined as the ratio of the 275 couplings for ti and bb, respectively, 
relative to SM-like strength. We have also taken into account the contributions of the 
"difficult" decay modes i) — viii) to the total decay width when kinematically allowed. 

The expected statistical significances for the processes 1) gg — > h/a — > 77 and 6) gg — > 
h -> WW« -> £ + t vv are computed from results for the SM Higgs boson taken from 
Ref. 1 8 1, Fig. 1. The application of the NMSSM corrections using Ri, ti and 6, (which 
determine T(gg — > hi), BR(hi — ► 77) and BR(hi — > WW*)) is straightforward in these 
two cases. The expected statistical significances for process 2) Wh/a or tih/a with r yy£ ± 
in the final state, are taken from the same figure. One finds that Whi and tihi production 
contribute with roughly equal weight to the SM signal. This allows us to decompose the 
expected significance into the corresponding production processes, apply the NMSSM 
corrections, and then recombine the production processes. 

For the Standard Model Higgs boson, the expected statistical significances for process 
3) tih/a with h/a — > bb, are taken from table 19-8 in Ref. 0, with the extension to Higgs 
boson masses above 120 GeV as provided by |15]. For the standard model process 5) 
gg -> h -> ZZ& -> 4 leptons, we again use Ref. |9[, tables 19-18 and 19-21. In both 
cases, the application of the NMSSM corrections is straightforward. The estimation of 
the statistical significances for the process 4) bbh/a with h/a — > t + t~, is more involved. 
Figure 19-62 of Ref. give the 5a contours in the tan (3 — tua plane of the MSSM. The 
critical issue is how much of these 5a signals derive from gg — ► H + gg — > A production 
and how much from associated bbH + bbA production, and how each of the gg fusion and 
bb associated production processes are divided up between H and A. For more details on 
our procedure, see Ref. |5|. 

Results for the statistical significances of the hi signals in modes 7) WW — > h — ► t + t~ 
and 8) WW — > h — > WW^*^ were similarly obtained by rescaling the theoretical results 

of ED. 

In the case of Higgs states close in mass, the individual statistical significances have to 
be recombined. We used the procedure of Ref. |16[, which gives the combined statistical 
significance N SD ^ 2 as a function of the individual statistical significances N S d,i, N S d,2 for 
two states with masses m 1? m 2 if \m\ — m 2 \ < Afc " 11 ^"" 12 , where A k is the mass resolution 
of the decay mode k): 

N 2 _ N 2 , N 2 o N S D,lN SD ,2 

' V A fe mi+m 2 / 
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For the processes 1) - 8), we assumed Ai )2) 5,6 = 1%, A 3i7i8 = 10%, and A 4 = 15%. 

Using the above procedures, for each point in the parameter space of the NMSSM we 
obtain the statistical significances predicted for an integrated luminosity of 100 fb" 1 for 
each of the detection modes 1) - 8). In order to obtain the statistical significances for the 
various detection modes at 300 fb" 1 , we multiply the 100 fb" 1 statistical significances by 
y/3 in the cases 1), 2), 3), 5) and 6), but only by a factor of 1.3 in the cases 4), 7) and 8). 
That such a factor is appropriate for mode 4), see, for example, Fig. 19-62 in |9[. Use of 
this same factor for modes 7) and 8) is simply a conservative guess. 

In our set of randomly scanned points defined in section |2j we selected those for 
which all the statistical significances [including VTl^-fusion modes 7) and 8)] are below 
5a. We obtained a lot of points, all with similar characteristics. Namely, in the Higgs 
spectrum, we always have a very SM-like CP-even Higgs boson with a mass between 115 
and 135 GeV (i.e. above the LEP limit), which can be either hi or h 2 , with a reduced cou- 
pling to the gauge bosons Ri ~ 1 or R 2 ~ 1, respectively. This state decays dominantly 
to a pair of (very) light CP-odd states, a\a\, with m ai between 5 and 65 GeV 1 . The singlet 
component of a\ has to be small in order to have a large hi — > aia\ or h 2 — > a\a\ branching 
ratio when the hi or h 2 , respectively, is the SM-like Higgs boson. Further, when the hi or 
h 2 is very SM-like, one has R' n ~ or R' 21 ~ 0, respectively, so that the e + e~ — > h-\_a\ 
or e + e~ — > h 2 ai associated production rate is very small and these processes thus place 
no constraint on the light CP-odd state at LEP. We have selected six difficult benchmark 
points. For points 1 - 3, hi is the SM-like CP-even state, while for points 4 - 6 it is h 2 . We 
have chosen the points so that h X 2 and a x have different masses. The main characteristics 
of the benchmark points are displayed in table [TJ Note the large branching ratios of the 
SM-like Higgs boson (hi for points 1-3 and h 2 for points 4 -6) to decay to axc^. Branching 
ratios for the non-SM-like h 2 (points 1 - 3) or hi (points 4-6) are also shown. For points 
4-6, with < 100 GeV, the hi is mainly singlet. As a result, R' n is very small and thus 
there are no LEP constraints on the hi and a x from e + e — > hid\ production. 

In the case of the points 1-3, the tabulated branching ratios of the h 2 are all very 
small. Even though the h 2 is heavy, it has very weak WW, ZZ coupling and BR(h 2 — > 
ZZ + WW) ~ 0.01. Further, the h 2 production rate due to gg fusion is greatly suppressed, 
as indicated in the table, as is its production rate from WW, ZZ fusion. This means that 
searches for the h 2 in the ZZ — > 4£ final state are impossible. In fact, the h 2 decays primar- 
ily to it with the 2nd most important mode being Zai. The only search channel for the 
h 2 that might be worth further investigation appears to be the bbh 2 Yukawa "radiation" 
process with h 2 — > it. Even though the bbh 2 Yukawa couplings are somewhat enhanced, 
the overall production rate will not be that large and we anticipate that the signal for the 
h 2 in this mode would be quite marginal. Regarding the LC, since m h2 >500 GeV, y/s e + e - 
substantially above the 500 GeV level would be required. Even for such y/s e + e -, h 2 de- 
tection would be very challenging since the ZZh 2 , WWh 2 couplings are very suppressed 
and the bbh 2 Yukawa coupling is only somewhat enhanced. Thus, we have focused on 
the aidi decays of the hi for points 1-3. 

For points 4-6, even though the h\ is quite light, it is dominantly singlet-like and 
therefore has suppressed couplings to all relevant SM particles, WW, ZZ, bb and it. As 
a result, it has very low production rates and its decays (primarily to either bb, points 

1 In some rare cases, one has hi decaying mainly to h\h\, but then the LHC significances corresponding 
to the discovery channels 1) — 8) are quite close to 5, so we have not considered cases of this type. 
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Point Number 


1 


2 


3 


4 


5 


6 


Bare Parameters 




A 


0.2872 


0.2124 


0.3373 


0.3340 


0.4744 


0.5212 


K 


0.5332 


0.5647 


0.5204 


0.0574 


0.0844 


0.0010 


tan/3 


2.5 


3.5 


5.5 


2.5 


2.5 


2.5 


A^cff (GeV) 


200 


200 


200 


200 


200 


200 


A x (GcV) 


100 





50 


500 


500 


500 


A K (GeV) 




















CP-even Higgs Boson Masses and Couplings 




m hl (GeV) 


115 


119 


123 


76 


85 


51 


Ri 


1.00 


1.00 


-1.00 


0.08 


0.10 


-0.25 


ti 


0.99 


1.00 


-1.00 


0.05 


0.06 


-0.29 


h 


1.06 


1.05 


-1.03 


0.27 


0.37 


0.01 


Relative gg Production Rate 


0.97 


0.99 


0.99 


0.00 


0.01 


0.08 


BR(hi -» bb) 


0.02 


0.01 


0.01 


0.91 


0.91 


0.00 


BR{hi -» r+r") 


0.00 


0.00 


0.00 


0.08 


0.08 


0.00 


BR(hi — > aiOi) 


0.98 


0.99 


0.98 


0.00 


0.00 


1.00 


m ft2 (GeV) 


516 


626 


594 


118 


124 


130 




-0.03 


-0.01 


0.01 


-1.00 


-0.99 


-0.97 


*2 


-0.43 


-0.30 


-0.10 


-0.99 


-0.99 


-0.95 


&2 


2.46 


-3.48 


3.44 


-1.03 


-1.00 


-1.07 


Relative gg Production Rate 


0.18 


0.09 


0.01 


0.98 


0.99 


0.90 


BR(h 2 -> &&) 


0.01 


0.04 


0.04 


0.02 


0.01 


0.00 


BR(h 2 — » t+t~) 


0.00 


0.01 


0.00 


0.00 


0.00 


0.00 


BR(h 2 — > aiOi) 


0.04 


0.02 


0.83 


0.97 


0.98 


0.96 


(GeV) 


745 


1064 


653 


553 


554 


535 


CP-odd Higgs Boson Masses and Couplings 




m Ql (GeV) 


56 


7 


35 


41 


59 


7 


*'i 


0.05 


0.03 


0.01 


-0.03 


-0.05 


-0.06 


&i 


0.29 


0.34 


0.44 


-0.20 


-0.29 


-0.39 


Relative gg Production Rate 


0.01 


0.03 


0.05 


0.01 


0.01 


0.05 


Bii(oi -> 66) 


0.92 


0.00 


0.93 


0.92 


0.92 


0.00 


Bii(oi -> r+r") 


0.08 


0.94 


0.07 


0.07 


0.08 


0.90 


m a2 (GeV) 


528 


639 


643 


560 


563 


547 


Charged Higgs Mass (GeV) 


528 


640 


643 


561 


559 


539 


Most Visible Process No. 


2 (h{) 


2(/ii) 


8(/ii) 


2(h 2 ) 


8(h 2 ) 


8(h 2 ) 


Significance at 300 fb _1 


0.48 


0.26 


0.55 


0.62 


0.53 


0.16 



Table 1: Properties of selected scenarios that could escape detection at the LHC. In the table, Ri, ti and bi 
are the ratios of the hi couplings to VV, ti and bb, respectively, as compared to those of a SM Higgs boson 
with the same mass; t[ and b[ denote the magnitude of the 175 couplings of ai to ti and bb normalized 
relative to the magnitude of the ti and bb SM Higgs couplings. We also give the production for gg — > hi 
fusion relative to the gg fusion rate for a SM Higgs boson with the same mass. Important absolute branching 
ratios are displayed. For points 2 and 6, BR(a\ — > jj) ~ 1 — BR(a\ — > t+t~). For the heavy h 3 and a 2 , we 
give only their masses. In the case of the points 2 and 6, decays of a\ into light quarks start to contribute. 
For all points 1-6, the statistical significances for the detection of any Higgs boson in any of the channels 
1) - 8) (as listed in the introduction) are tiny; their maximum is indicated in the last row, together with the 
process number and the corresponding Higgs state. 
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4 and 5, or a^ax, point 6) will not produce particularly distinctive final states, especially 
given the smallness of m hl . Thus, we believe the hi would be very difficult to detect at 
the LHC. At the LC, the hi has ZZ coupling sufficient for detection only in the case of 
point 6. Given the complicated nature of the hi — ► aiai final state in this case (including 
the fact that the ai is so light that ai — ► r + r is dominant) detection would have to rely 
on the ZX reconstructed recoil mass technique. Thus, for points 4-6, we will focus on 
searching for the SM-like h 2 in its dominant h 2 — > a x ai decay mode. 

Of course, we should also consider whether or not the h 3 or a 2 might be detectable 
for points 1-6. Since m h3 is very large in all cases, one would need to rely on the h 3 — > 
ZZ — > A£ mode [mode 5) in our first list] for its detection. But, the gg and WW fusion 
production rates are suppressed and BR(h 3 — > ZZ) < 0.01. Mode 5) for the h 3 contributes 
negligibly to the net statistical significance in the last column of table HI At the LC, ^/s e + e - 
substantially above 500 GeV would be required at all points. Detection of the a 2 would 
appear to be even more difficult. If a super high energy LC is eventually built, the pair 
production modes such as e + e — ► Z — > h s a 2 will be substantial and produce very viable 
signals. 

In the case of points 2 and 6, it should be noted that the a x — > t + t decays are domi- 
nant, with ai — > jj decays making up most of the rest. For points 1 and 3-5 for which 
BR(ai — > bb) is substantial, the b jets will not be that energetic and tagging will be some- 
what inefficient. Thus, we have chosen not to implement 6-tagging as part of the experi- 
mental procedures detailed in the next section. 

Finally, we should consider whether or not a direct search for the ai might be feasible. 
For our points 1-6, the a\ has suppressed couplings to both bb and tb and will thus be very 
weakly produced via Yukawa radiation and would not be detectable in this way at either 
the LHC or the LC IT7I . The only processes that would have significant rate are those 
relying on the ZZa\(i\ and WWa\d\ quartic couplings required by gauge invariance. At 
the LHC, the rate for, say, WW — > a\a x is not exactly large and, as we shall see, there 
is a large background from it production. Rates for continuum WW — > oiai production 
at the LC are substantial for a light a\ |18J, but again we shall see that there is a large 
background. In the next section, we demonstrate that the constraint that M aiai ~ m h with 
mh ~ 100 GeV is absolutely critical in order to have any hope of extracting a signal in 
WW — > h production at both the LHC and the LC. 

4 New Channels at the LHC and LC 

As we have already stressed, for the points summarized in table [T] the a\ is light and de- 
cays almost entirely into bb (or jj for points 2 and 6) and t + t~. The possible final states 
are thus bbbb (or Aj for points 2 and 6), bbr + T~ (or 2jr + r~) and t + t~t + t~ . A 46-signal 
would be burdened by a large QCD background even after implementing 6-tagging. 
A 4j-signal would be completely swamped by QCD background. Meanwhile, the 4r- 
channel would not allow one to reconstruct the hi, h 2 resonances. Hence, in attempting 
to cover these problematic points of the NMSSM parameter space at the LHC, we will 
focus in this study on the 2b2r (or 2j2r) signature. The extra factor of two in BR(aiai — > 
bbr+T-) = 2BR(ai -> bb)BR(a 1 -> t + t~) relative to BR(a iai -> bbbb) = [BR(a 1 -»■ bb)} 2 
and BR(ai<ii — ► t + t~t + t~) = [BR(ai — ► t + t~)] 2 (and similarly with bb — ► jj for points 
2 and 6) is an additional advantage. It compensates in part (or further increases) the 
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Yukawa suppression (or enhancement) with respect to the hi/h 2 — > aidi — > 46(4r) decay 
rate. 2 In addition, we will be looking at r's decaying leptonically to electrons and muons, 
yielding some amount of missing (transverse) momentum, p^ iss , that could be projected 
onto the visible e, /^-momenta in an attempt to reconstruct the parent r-direction. Since 
for points 2 and 6 the a\ does not decay to bb and since the b and b that do come from 
ai are not very energetic given the modest m ai mass for points 1 and 3-5, we will not 
employ 6-tagging as part of our analysis. 

4.1 Results for the LHC 

Detection of WW* and t + t~ decays of a relatively light Higgs boson with SM-like WW 
coupling is straightforward at the LHC in the vector- vector fusion production mode [i.e. 
the Higgs production channel entering the signatures 8) and 9)] when the light Higgs bo- 
son has mass above the LEP limits and decays directly to t + t (the status of the bb mode 
has not yet been determined). Because the significance of this type of signal for a SM-like 
Higgs boson is very large, we consider here the same W^VT-fusion production mode in the 
context of the NMSSM. (We reemphasize that the hi [cases 1 - 3] or h 2 [cases 4-6] has 
nearly full SM strength coupling to WW.) However, the bbr + T~ (or 2jt + t~, for points 
2 and 6) final state of relevance is more complex and subject to larger backgrounds than 
is a simple t + t final state. Further, the a\ masses of interest are only half as large as 
the Higgs masses for which the direct h — ► t + t - signals are viable. In order to extract 
the 2j2r NMSSM Higgs boson signature from the central detector region, we have ex- 
ploited forward and backward jet tagging on the light quarks emerging after the double 
VT-strahlung preceding VTiy-fusion, the utility of which is reviewed in Ref. lli"9"ll . If we 
require two additional forward/backward jets, it is clear that the leading background is 
due to tt production (since we are assuming a heavy SUSY spectrum) and decay via the 
purely SM process, 

gg ^ tt - bbW + W~ - bbr+T- + pl iss , (8) 

in association with forward and backward jet radiation. 
In summary, at the LHC, the signature is 

• 2 forward/backward jets, at least 2 central jets, p^ iss and a t + t~ pair decaying lep- 
tonically (to electrons and/ or muons). 

In order to carry out realistic numerical simulations, we have used a modification of the 
MSSM implementation [20[ of the HERWIG event generator |21|, in conjunction with the 
GET JET code (22] for calorimeter emulation and jet reconstruction. An ISAWIG format 
[23 1 input file has been edited by hand to incorporate the Higgs boson mass spectrum 
and decay rates as predicted for each of the NMSSM points 1-6, while in the main 
HERWIG code (v6.4) the subroutine implementing the vector-vector fusion process has 
been modified to account for the different Higgs- vertices pertaining to the NMSSM 
points 1-6 given in table HJ The above codes do not include K factors for either the signal 
or the background. 

An outline of the selection procedure and cuts used is the following. 

2 Recall that it is the running 6-quark mass (evaluated at the a\ mass) that enters the Yukawa coupling, 
rather than the pole mass. 
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• Acceptance cuts: 

l^jetl < 5, pf et > 20 GeV, Ai? jet _ jet > 0.7, 

^r^S n <0, he P ton\ < 2.5, 

Pkpton > 10 GeV, no lepton isolation. 

• Since the ai will not have been detected previously, we must assume a value for 
m ai . It will be necessary to repeat the analysis for densely spaced m ai values and 
look for the m ai choice that produces the best signal. 

• We look among the central jets for the combination with invariant mass Mjj closest 
to m ai (no 6-tagging is enforced, b's are identified as non-forward/backward jets). 

• Select the two highest transverse-momentum leptons in any flavor combination and 
with opposite charge. After ensuring that these are not back-to-back (by requiring 
that their relative angle is smaller than 175 degrees), resolve the p^ iss along their 
directions and reconstruct the invariant mass, M T + T - . 

• Plot the Mjj T + T - invariant mass using the four four-momenta reconstructed in the 
two previous steps, as seen in the top plot of Fig. [T] — the plot presented assumes 
that we have hit on the correct m ai choice. 

The selection strategy adopted is clearly efficient in reconstructing the hi (for points 1- 
3) and h 2 (for points 4-6) masses from the jjr + r~ system, as one can appreciate by noting 
the peaks appearing in the LHC plot of Fig.[l]at M jjT + T - m 100 GeV. In contrast, the heavy 
Higgs resonances at mh 2 for points 1-3 and the rather light resonances at for points 
4-6 (recall table [TJ> do not appear, the former mainly because of the very poor production 
rates and the latter due to the fact that either the hi — > a\ai decay mode is not open (points 
4, 5) or - if it is - the 6-quarks and e/yu-leptons eventually emerging from the ai decays are 
too soft to pass the acceptance cuts (point 6, for which m ai = 7 GeV and m hl =51 GeV). 
For all six NMSSM setups, one may see a hint of a resonance in the data in the very end of 
the low mass tail of the tt background (see the insert in the top frame of Fig. [TJ. However, 
after summing the background distribution and any one of the signal spectra, it could 
be difficult to ascertain the existence of the h\ or h 2 peaks from the net line shapes. Still, 
statistics are significant. To estimate S/y/B, we assume L = 300 fb" 1 , a K factor of 1.1 for 
WW fusion and a K factor of 1.6 for the it background. (These K factors are not included 
in the plots of Fig. [T]) We sum events over the region 60 < Mj jT + T - < 90 GeV. For points 
1-6 we obtain signal rates of about S = 890, 600, 750, 1030, 915, 500, respectively. The tt 
background rate is B ~ 320. This gives N SD = S/VB of 50, 34, 42, 58, 51, 28 for points 
1-6, respectively. These are substantial. However, given the broad distribution of the 
signal, it is clear that the crucial question will be the accuracy with which the background 
shape can be predicted from theory. The background normalization after the cuts imposed 
in our analysis would be very well known from the higher M JJT + T regions. 

4.2 The LC scenario 

While further examination of and refinements in the LHC analysis may ultimately lead 
us to have full confidence in the viability of the NMSSM Higgs boson signals discussed 
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Figure 1: Reconstructed mass of the jjr + r system for signals and backgrounds after the selections 
described, at the LHC (top) and a LC (bottom). We plot da/dM jjT+T - [fb/10 GeV] vs M jjT+T - [GeV] using 
GET JET and HERWIG 6 . 4 with: IPROC=372 adapted to NMSSM couplings and decay rates for the signal 
and IPROC=1706 for the background at the LHC; IPROC=920 adapted to NMSSM couplings and decay 
rates for the signal and IPR0C=12 6 (250) for the tt[ZZ] background at a LC). Statistics used: 500,000 
points. Normalization is to the total cross section after cuts. In both plots, the lines corresponding to points 
4 and 5 are visually indistinguishable. No K factors are included. 
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above, an enhancement at low M JJT + r of the type shown (for some choice of m ai ) will 
nonetheless be the only evidence on which a claim of LHC observation of Higgs bosons 
can be based. Ultimately, a means of confirmation and further study will be critical. Thus, 
it is important to summarize the prospects at the LC, with energy up to 800 GeV, for 
detecting the NMSSM Higgs bosons, especially in the context of the difficult scenarios 1 — 
6 of table[T]discussed here. As we summarize below, the LC is certainly guaranteed to find 
the hi (points 1-3) or the h 2 (points 4-6) and will allow a good determination of the h%aiai 
or h 2 aiai branching ratio, respectively, a possibly very important piece of information 
for unraveling the model. In these scenarios, the LC would be an absolutely essential 
complement to the LHC, and, in particular, would be the only machine at which precision 
Higgs physics studies could be pursued. In what follows, we will use the notation h for 
h = h\ for points 1-3 and h = h 2 for points 4-6 in table [TJ 

Because the ZZh coupling is nearly full strength in all cases, and because the h mass 
is of order 100 GeV, discovery of the h will be very straightforward via e + e~ — > Zh using 
the e + e — ► ZX reconstructed Mx technique which is independent of the "unexpected" 
complexity of the h decay to aiai. This will immediately provide a direct measurement 
of the ZZh coupling with very small error |24|. The next stage will be to look at rates 
for the various h decay final states, F, and extract BR{h — > F) — a(e + e~ — > Zh — ► 
ZF)/a(e + e~ — > Zh). For the NMSSM points considered here, the main channels would 
be F = bbbb, F = bb~T + T~ and F = t + t~t + t~ ? At the LC, a fairly accurate determination 
of BR{h — > F) should be possible in all three cases. This is important if we are to be able 
to determine BR{h — > a x a x ) independently. The procedure is the following. From 

BR(h -> bbbb) = BR(h -> a x a x )[BR{a x -> bb)] 2 
BR(h^T + T-T + r-) = BR(h^ ai ai)[BR(ai^T + r-)] 2 

BR(h — > b~b~T + T~) = 2BR(h — > aiai )BR(a x -> bb)BR(ai -> t + t~) , (9) 

we see that BR(ai — > t + t~) / BR(a>i — > bb) can be determined using two independent 
ratios. Assuming that these are the only two decay channels (as can be checked by looking 
for other final states), the requirement BR(ai — > t + t~) + BR(ai — > bb) = 1 then allows 
extraction of the individual BR's. Once these are known, we can use the above equations 
to determine BR(h — > aiai). Certainly, this procedure will be carried out in the e + e~ — > 
Z/i production mode. However, we have not performed the required simulation for this 
paper to determine the actual accuracy that can be achieved. 

Instead, we have considered the equally (or perhaps more) useful vector-vector fusion 
mode that will be active at a LC. In fact, at 800 GeV or above, it is the dominant Higgs 
boson production channel for CP-even Higgs bosons in the intermediate mass range. 
Contrary to the case of the LHC though, the dominant contribution (from WW fusion) 
does not allow for forward and backward particle tagging, as the incoming electron and 
positron convert into (anti)neutrinos, which escape detection. Although the ZZ fusion 
contribution would allow tagging of forward/backward e~ and e + , the cross section is a 
factor of 10 smaller (see Fig. 4 of Ref. |24|) in comparison. Still, since the it background 
would not typically be accompanied by forward /backward e + and e~ in e + e~ collisions, 
unless initiated by 77 interactions (which tend to produce less energetic it pairs), tagging 
of the electron and positron in the final state might offer a handle against the it noise. 

3 Here, and in the following discussion, bb should be replaced by jj for points 2 and 6, where j refers to 
any possible non-6 jet. 
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But, in any case it seems unnecessary to focus on ZZ fusion given that it production at 
the LC proceeds through EW interactions, rather than via QCD as at the LHC. Indeed, in 
going from the LHC to a LC, the it production rate should be reduced (apart from parton 
distribution functions) by a factor of 0(a s /a cm ) 2 with respect to the vector boson fusion 
Higgs boson signal (which is an electroweak process at both the LC and the LHC). A new 
background of importance does emerge at the LC, namely ZZ-pair production with one 
Z decaying to jj and the other to t + t pairs. The ZZ background is negligible at the 
LHC mainly because associated energetic forward /backward jets from QCD initial state 
radiation are infrequent. At the LC, the ZZ background plays a more significant role and 
has been simulated in our HERWIG and (LC-adjusted) GET JET numerical analysis. 

At a LC, the optimal signature will thus be different than at the LHC and a different 
set of selection criteria are needed. We choose selection criteria that will retain both the 
WW and the ZZ fusion Higgs boson production processes. 

• We require an even number of oppositely-charged leptons (n t = 2, 4...). 

• If n>£ > 4 then we demand that two of these must be an electron and a positron, so 
that the final state would be consistent with having been generated by forward /back- 
ward e ± tagging in ZZ-fusion and t + t decays in which both r's decay leptonically 
(to electron and muons). 

• Finally, we require at least 2 central jets and significant p^ iss . 

Given the required final state as above, we employ the following selection procedure and 
cuts. 

• Acceptance cuts: 

| cos^ e t | < 0.990, pf ct > 5 GeV, Ai? jct _ jet > 0.4, 

C+ x • V?- n < (if n t > 4), | cos # lepton | < 0.995, 
Pkpton > 5 GeV, no lepton isolation. 

• We look among the central jets for the combination with invariant mass Mjj closest 
to m ai (again, no 6-tagging is enforced — b's are identified as non-forward/backward 
jets). 

• Out of the n e leptons, upon excluding the e + e~ pair in which the electron and 
positron are those with the largest rapidities if n e > 4, select the two with high- 
est transverse-momenta in any flavor combination and with opposite charge. 

After ensuring that these are not back-to-back, resolve the p^ liss along their directions 
and reconstruct the invariant mass M T + T -. 

• Plot the Mjj T + T - invariant mass (see the bottom of Fig.[T} using the four four-momenta 
reconstructed in the two previous steps. 

Note that it is not fruitful to place cuts on the invariant masses Mjj and M T + T - that ex- 
clude Mjj, M T + T - ~ mz in an attempt to reduce the ZZ background. This is because 
the SM-like h mass is typically of order 115 GeV, i.e. not so far from mz, and the experi- 
mental resolutions in the two masses Mjj and M T + T - are poor, either because of the large 
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number of hadronic tracks or the missing longitudinal momenta of the (anti) neutrinos, 
respectively. 

Just as at the LHC, the peaks corresponding to the Higgs boson with weak VV cou- 
pling (h 2 for points 1-3 and hi for points 4 - 6) do not appear in the M jjT + T - distribution 
which can be reconstructed at a LC. Here, though, for the visible h (hi for points 1-3, 
h 2 for points 4-6) the situation is better than at the LHC. For all points 1-6 the mass 
peaks (again centered at 100 GeV) are now very clearly visible above both the it and ZZ 
backgrounds, particularly for the case of points 2 and 6 (see insert in the bottom frame 
of Fig. [!}. Assuming L = 500 fb _1 , the points 1,3,4,5 yield 5 events per 10 GeV bin, on 
average, in the 50 to 150 GeV mass interval of interest. This would constitute a convinc- 
ing signal given the very small size predicted for the background. Notice that, although 
to assign the entire missing transverse momentum to the r-lepton system may seem not 
entirely appropriate (as it is largely due to forward/backward (anti)neutrinos from the 
incoming electrons and positrons in WW fusion), this does not hamper the ability to re- 
construct the Higgs mass peaks. However, this implies that a proportion of the signal 
events will tend to reproduce the overall yjs e + e - value in the Mj jT + T - distribution. The 
effect is more pronounced for points 1 and 3-5, which is where the a x mass is larger 
(see Table so that most of the hadronic tracks composing the emerging jets easily en- 
ter the detector region. For points 2 and 6, where m ai is below 10 GeV, this may often 
not be true and it appears that the consequent effect of these hadrons escaping detec- 
tion is that of counterbalancing the p^ iss contributions related to the neutrinos left behind 
in WW fusion reactions. For the case of the ZZ noise, in the presence of full coverage 
and perfect resolution of the detector, one would have Mjj T + T - = Js e + e -, which explains 
the very noticeable preference for this process to produce a concentration of events with 
Mj jT + T - around 800 GeV in the mass distribution. (The "tails" beyond ^Js e + e - are due to 
the smearing of the visible tracks in our Monte Carlo analysis.) 

Finally notice that we have included Initial State Radiation (ISR) and beam-strahlung 
effects, as predicted using the HERWIG default. These tend to introduce an additional 
unresolvable missing longitudinal momentum, although to a much smaller extent than 
do the Parton Distribution Functions (PDFs) in hadron-hadron scattering at the LHC. 

A final note regarding WW, ZZ fusion to h. Since we will know from Zh production 
the magnitude of the WWh and ZZh coupling, the rate for e + e~ — > ZZ, WW — > h — > 
bbr + T~ will determine BR[h — > 66r + r~), with similar results perhaps possible for other 
channels. Thus, we could hope to implement the procedure described earlier for extract- 
ing BR(h — > ai<2i) using both WW, ZZ — > h fusion and Zh production. 

5 Conclusions 

In summary, if the NMSSM parameters are such that the most SM-like of the CP-even 
Higgs bosons, h, is relatively light and decays primarily to a pair of CP-odd Higgs states, 
h — > aa, then there will be a statistically highly significant LHC signal (from WW — ► h — > 
aa) of an S/B ~ (500 - 1000) /300 bump (for L = 300 ftT 1 ) in the low-mass tail of a rapidly 
falling jjr + r~ mass distribution. The LHC would thus give a very strong indication of the 
presence of a Higgs boson. However, this detection mode is not exactly in the gold-plated 
category. The LC will be absolutely essential in order to confirm that the enhancement 
seen at the LHC really does correspond to a Higgs boson. At the LC, discovery of a light 
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SM-like h is guaranteed to be possible in the Zh final state using the recoil mass technique. 
Further, we have seen that WW, ZZ fusion production of the h will also produce a viable 
signal in the jjr + r~ final state (and perhaps in the Aj and t + t~t + t~ final states as well, 
although we have not examined these |25[). 4 Measurement of the relative rates for the 
Aj, 2j2r and At final states will allow a determination of the a — > jj and a — > t + t~ 
branching ratios, which in turn will make possible the determination of BR{h — > act), 
a potentially very important measurement. Unfortunately, the standard techniques for 
determining the total width of a SM-like h relying on the 77 and/ or WW final state decays 
of the h will not be available and it is thus unlikely that we can convert a measurement 
of BR{h — > aa) into a determination of the partial width F(h — > aa) (which would be 
most directly related to the very interesting haa coupling strength). As we have stressed, 
for parameter space points of the type we have discussed here, detection of any of the 
other MSSM Higgs bosons is likely to be impossible at the LHC and is likely to require 
an LC with ^s e + e - above the relevant thresholds for h'a' production, where hi and a! are 
heavy CP-even and CP-odd Higgs bosons, respectively. Although results for the LHC 
indicate that Higgs boson discovery will be challenging for the type of situations we have 
considered, improved techniques for extracting a signal are likely to be developed once 
data is in hand and the ti background can be more completely modeled. Clearly, if SUSY 
is discovered and no Higgs bosons are detected in the standard MSSM modes, a careful 
search for the signal we have considered should have a high priority. Finally, we should 
remark that the h — > aa search channel considered here in the NMSSM framework is also 
highly relevant for a general two-Higgs-doublet model, 2HDM. It is really quite possible 
that the most SM-like CP-even Higgs boson of a 2HDM will decay primarily to two CP- 
odd states. This is possible even if the CP-even state is quite heavy, unlike the NMSSM 
cases considered here. 
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